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Abstract 

Order positions are key variables in algorithmic trading. This paper studies the limiting 
behavior of order positions and related queues in a limit order book. In addition to the fluid and 
diffusion limits for the processes, fluctuations of order positions and related queues around their 
fluid limits are analyzed. As a corollary, explicit analytical expressions for various quantities of 
interests in a limit order book are derived. 


1 Introduction 

In modern financial markets, automatic and electronic order-driven trading platforms have largely 
replaced the traditional floor-based trading; orders arrive at the exchange and wait in the Limit 
Order Book (LOB) to be executed. There are two types of buy/sell orders for market participants 
to post, namely, market orders and limit orders. A limit order is an order to trade a certain amount 
of security (stocks, futures, etc.) at a given specified price. Limit orders are collected and posted 
in the LOB, which contains the quantities and the price at each price level for all limit buy and sell 
orders. A market order is an order to buy/sell a certain amount of the equity at the best available 
price in the LOB; it is then matched with the best available price and a trade occurs immediately 
and the LOB is updated accordingly. A limit order stays in the LOB until it is executed against a 
market order or until it is canceled; cancellation is allowed at any time without penalty. 

The availability of both market orders and limit orders presents market participants opportu¬ 
nities to manage and balance risk and profit. As a result, one of the most rapidly growing research 
areas in financial mathematics has been centered around modeling LOB dynamics and/or minimiz¬ 
ing the inventory/execution risk with consideration of the microstructure of LOB. A few examples 
include laiaiSEiiiTiiisiEainiEsiETiisiiMiiiiiiiziiiaiiTiEoj. 

At the core of these various optimization problems is the trade-off between the inventory risk 
from unexecuted limit orders and the cost from market orders. While it is straightforward to 
calculate the costs and fees of market orders, it is much harder to assess the inventory risk from 
limit orders. Critical to the analysis is the dynamics of an order position in an LOB. Because of the 
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price-time priority (i.e., best-priced order first and first-in-first-out) in most exchanges in accordance 
with regulatory guidelines, a better order position means less waiting time and a higher probability 
of the order being executed. In practice, reducing low latency in trading and obtaining good order 
positions is one of the driving forces behind the technological race among high-frequency trading 
firms. Recent empirical studies by Moallemi and Yuan [H] show that values of order positions (if 
appropriately defined) have the same order of magnitude of a half spread. Indeed, analyzing order 
positions is one of the key components for studying algorithmic trading strategies. Knowing both 
the order position and the related queue lengths not only provides valuable insights into the trading 
direction for the “immediate” future but also provides additional risk assessment for the order — 
if it were good to be in the front of any queue, then it would be even better to be in the front of a 
long queue. Therefore, it is important to understand and analyze the dynamics of order positions 
together with their related queues. This is the focus of our work. 


Queue Tail 
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r" 
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Figure 1: Orders happened in the best bid queue. 


Our contributions. The dynamics of an order position in a queue will be affected by both the 
market orders and the cancellations, and the dynamics of its relative position in a queue will be 
affected by the limit orders as well (see Figure 1). Without loss of generality, we will focus on an 
order position in the best bid queue along with the best bid and ask queues. Order positions in 
other queues will be similar and simpler because of the absence of market orders. 

First, we derive the fluid limit for the order positions and related best bid and ask queues; in a 
sense, this is a first order approximation to the processes. We show (Theorem IIII and Theorem [31]) 
that the rate of the order position approaching zero is proportional to the mean of order arrival 
intensities and to the summation of the average size of market orders and the “modified” average size 
of cancellation orders in the queue; this modification depends on different assumptions on order 
cancellations. We also derive the (average) time it takes for the order position to be executed. 
The derivation is via two steps. The first step is to establish the functional strong law of large 


2 




numbers for the related bid/ask queues; this is straightforward. The second step is intuitive but 
requires a delicate analysis involving passing the convergence relation of stochastic processes in 
their corresponding cadlag space with the Skorokhod topology to their integral equations. 

Next, we proceed to the second order approximation for order positions and related queues. 
The first step is to establish appropriate forms of the diffusion limit for the bid and ask queues. We 
establish a multi-variate functional central limit theorem (FCLT) using ideas from random fields. 
Under appropriate technical conditions, we show fTheorem ll4l) that the queues are two-dimensional 
Brownian motion with mean and covariance structure explicitly given in terms of the statistics of 
order sizes and order arrival intensities. The second step is to combine the FCLTs and the fluid 
limit results to show (Theorem 1 15 jl that fluctuations of the order positions are Gaussian processes 
with “mean-reversion”. The mean-reverting level is essentially the fluid limit of order position 
relative to the queue length modified by the order book net flow, which is defined as the limit order 
minus the market order and the cancellation. The speed of the mean-reversion is proportional to 
the order arrival intensity and the rate of cancellations. 

Our results are built on fairly general technical assumptions (stationarity and ergodicity) on 
order arrival processes and order sizes. For instance, order arrival processes (Section S]) can be 
Poisson processes or Hawkes processes, both of which have been extensively used in LOB modelings; 
see for instance, Abergel and Jedidi [2] and Huang, Lehalle, and Rosenbaum [33] . 

Practically speaking, studying order positions gives more direct estimates for the “value” of 
order positions, which is useful in algorithmic trading. Indeed, based on the fluid limit, we derive 
(Section Hj) explicit analytical comparisons between the average time an order is executed and the 
average time any related queue is depleted. This is an important piece of information especially 
when combined with an estimate on the probability of a price increase. The latter is a core quantity 
for the LOB and has been studied in Avellaneda and Stoikov [6| and Cont and de Larrard |2n[ll9| for 
special cases. In addition, we derive from the fluctuation analysis explicit expressions for the first 
hitting times of the queue depletion, for the expected order execution time, and for the fluctuations 
of order execution time and first hitting times. Furthermore, by the large deviations theory, we 
derive the tail probability that the queues deviate from their fluid limits. 

Related work. The main idea behind our analysis is to draw connections between LOBs and 
multi-class priority queues, as LOBs with cancellations are reminiscent of reneging queues; see 
for instance Ward and Glynn [Ml Eg. In the mathematical finance literature, there have been a 
number of papers on modeling LOB dynamics in a queuing framework and establishing appropriate 
diffusion and fluid limits for queue lengths or the order book prices. This line of work can be traced 
back to Kruk [39], who established diffusion and fluid limits for prices in an auction setting and 
showed that the best bid and ask queues converge to reflected two-dimensional Brownian motion 
in the first quadrant. Similar results were later obtained by Gont and de Larrard m for the best 
bid and best ask queues under heavy traffic conditions, where they also established the diffusion 
limit for the price dynamics under the same “reduced form” approach with stationary conditions 
on the queue lengths [20] . Abergel and Jedidi [T] modeled the volume of the order book by a 
continuous-time Markov chain with independent Poisson order flow processes and showed that mid 
price has a diffusion limit and that the order book is ergodic. Horst and Paulsen [32] studied the 
fluid limit for the whole limit order books including both prices and volumes, under a very general 
mathematical setting. Their analysis was further extended in Horst and Kreher m, where the 
order dynamics could depend on the state of the LOB. Under different time and space scalings. 
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Blanchet and Chen [9] derived a pure jump limit for the price-per-trade process and a jump diffusion 
limit for the price-spread process. 

One of our results, Theorem [T^ is mostly related to yet different from the diffusion limit in 
m- This is a result of a different scaling approach. In order for us to analyze the dynamics of the 
order positions, we need to differentiate limit orders from market and cancellation orders, whereas 
in m order processes are aggregated from limit, market, and cancellations orders. Because of this 
aggregation, they could use the main idea from “heavy-traffic-limit” in classical queuing theory 
and assume that the mean order flow is dominated by the variance. While this assumption m 
Assumption 3.2] is critical to their analysis, it does not hold in our setting where each individual 
order type is considered. On the other hand, if we were to impose this assumption, then our result 
will be reduced to theirs because the second term in Eqn. (j3.7h would simply vanish. 

To the best of our knowledge, the dynamics of order positions and its relation to the queue 
lengths, which is the focus of our work, has not been studied before. Indeed, classical queuing 
tends to focus more on the stability of the entire system, rather than analyzing individual requests. 
Most of the existing modeling approaches in algorithmic trading have ignored order positions, with 
very limited efforts on the probability of it being executed. For instance, such a probability is 
either assumed to be a constant as in Cont and de Larrard [IHl [20] and Guo, de Larrard and 
Ruan [29], or is computed numerically from modeling the whole LOB as a Markov chain as in Hult 
and Kiessling [34], or is analyzed with a homogeneous Poisson process for order arrivals and with 
constant order sizes as in Cont, Stoikov, and Talreja [22]. 


2 Fluid limits of order positions and related queues 

2.1 Notation 


Without loss of generality, consider the best bid and ask queues. Then there are six types of orders: 
best bid orders (bb), market orders at the best bid (mbb), cancellation at the best bid (ebb), best 
ask (ba), market orders at the best ask (mba), and cancellation at the best ask (cba). Denote the 
order arrival process by N = > 0) with the inter-arrival times {Dj}i>i. Here 

{ m 

m : Di <t 

i=i 



For simplicity, assume that there are no simultaneous arrivals of different types of orders. Consider 
order arrivals of any of these six types as a point process, and define a sequence of six-dimensional 
random vectors {V where for the ith order 


:= (Vl, • • •, v!^), 


represents the sizes of the six types of orders; by the assumption, exactly one entry of 1/ j is positive. 
For instance, I^s = (0,0,0,4,0,0) means the fifth order is of size 4 and of type ba, i.e., a limit order 
at the best ask. In this paper, we only consider cadlag processes. 

For ease of references in the main text, we will use the following notation. 


• D[0,T] is the space of one-dimensional cadlag functions on [0,T], while D^[0,T] is the space 
of AT-dimensional cadlag functions on [0,T]. Consequently, the convergence in this space is, 
unless otherwise specified, in the sense of the weak convergence in D^[0, T] equipped with Ji 
topology; 
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• Loo[0,T] is the space of functions / : [0,T] —equipped with the topology of uniform 
convergence; 

• ACq[^,T] is the space of functions / : [0,T] —>■ that are absolutely continuous and 

/( 0 ) = 0 ; 

• ACq[^,T] is the space of non-decreasing functions / : [0, T] —)■ that are absolutely con¬ 

tinuous and /(O) = 0. 

Similarly, we define oo), D^[0,oo), Loo[0, oo), ^Co[0, oo), ^Cq"[0, oo) for T = oo. 


2.2 Technical assumptions and preliminaries 

In order to study the fluid limit for the order position and related queries, we will first need to 
impose some technical assumptions. 

Assumption 1. {Di}i>i is a stationary array of positive random variables with 

Di + D2 + --- + Di 1 , uhi-i 

- m probability 

i A 

as i ^ oo, where X is a positive constant. 

Assumption 2. i}i>i is a stationary array of square-integrable random vectors with 


V, 


in probability 




ii>l 


as i ^ oo, where V = iV^, V^,..., V^) is a constant vector. 
Assumption 3. is independent o/{I^j}j 

Now, we define a new process as follows, 

n 


Cn{t) = - V 7^. 

71 


2=1 


We call such a process the scaled net order flow process. 

Theorem 4. Given Assumptions\^ and\^ for any T > 0, 

'^n^XVe, in {D^[0,T], Ji) as n ^ oo, 
where e is the identity function. 

Proof. First, we define the scaled processes and by 

[nt\ 

n 

2 = 1 

\nt\ 




( 2 . 1 ) 


2 = 1 
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Then by Assumption [T] and according to Glynn and Whitt [251 Theorem 5], the strong Law of 
Large Numbers (SLLN) also follows, i.e., 


lim 

i^oo 


Di + D 2 + • • • + Di 1 


a.s. 


i A’ 

Then by the equivalence of SLLN and FSLLN [25l Theorem 4], it is clear that for any T > 0, 

1 e 

Di ^ , a.s. in (iA[0, T], Ji) as n ^ 00 . 

n A 


2=1 


Moreover, since is square-integrable, it follows that IE[1L/] < 00 for 1 < j < 6. Note that 
{V^'\i>i is stationary and applying Birkhoff’s Ergodic Theorem m Theorem 6.28] leads to 


1 "■ 

-EL^e|f/|pi. 


a.s. as n ^ 00 , 


i=l 


where is the invariant u-algebra of {V^}i>i. Given the WLLN for it follows that 

E[e/ I D] = 


and 

1 " 

— > V-^ ^ a.s. as n —)• 00 . 

n ^ * 

2=1 

Therefore, again by [251 Theorem 4], 

-> • 1 

14-^ ^ V^e, a.s. in (L*[0, T], Ji) as n —>■ 00 

^ i=i 

Since the limit processes for {S^}„>i and {Sp}n>i, 1 < j < 6, are deterministic, then according 
to [531 Theorem 11.4.5], 

S^) ^ ^Ee, , a.s. in (iA’^[0,T], Ji) as n ^ 00 . 

Finally, from [53l Theorem 9.3.4], 

^n=>AlAe, in (L)®[0, T], Ji) as n —)• 00 . □ 

Next we proceed to study the order position in the best bid and related queues of the best bid 
and the best ask. We further assume 

Assumption 5. Cancellations are uniformly distributed on every queue. 

We will see that this assumption on cancellation is not critical, except for affecting the exact 
form of the fluid limit for the order position. (See Theorem [31] without this assumption in Section 

13) 
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Now define the scaled queue lengths with for the best bid queue and Q“ for the best ask 
queue, and the scaled order position Z„ by 

' Qiit) = Qim+c^it) - ciit) - 

, Ql{t) = Q“(0) + - Cl{t) - Cl{t), ^2.2) 

dZn{t) = -dCl{t) - ^^^dClit). 

The above equations are straightforward: bid/ask queue lengths increase with limit orders and 
decrease with market orders and cancellations according to their corresponding order flow processes; 
an order position will decrease and move towards zero with arrivals of cancellations and market 
orders; new limit orders arrivals will not change this particular order position; however, the arrival 
of limit orders may change the speed of the order position approaching zero following Assumption 
[5l hence the factor of 

Strictly speaking, Eqn. (12.211 only describes the dynamics of the triple (Q^(t), <5“ (t), Z„(t)) 
before any of them hits zero: hitting zero means that the order placed has been executed, 

while Q“ hitting zero means that the best ask queues is depleted. Since our primary interest is in 
the order position, with little risk we may truncate the processes to avoid unnecessary technical 
difficulties on the boundary. That is, define 

Tn = min{<,r“,r^}, (2.3) 


with 


Tn = inf{i > 0 : < 0}, t“ = inf{t > 0 : Qf{t) < 0}, rf = inf{t > 0 : < 0}. 

Now, define the truncated processes 

Quit) = Quit Q“(t) = Arn), Zn{t) = Zn{t A Tn) ■ (2.4) 


Still, it is not immediately clear that these truncated processes would be well defined either since 
we do not know a priori if the term is bounded when hits zero. This, however, turns 

out not to be an issue. 


Lemma 6. Eqn. (12.4p is well defined, with Zn{t) < Qnit) for any time t < min(r^,r“). In 
particular, rf < r^. 


Proof. Note that is a positive jumping process. Therefore, when 5Zn{t) = 0, we have 6C^{t) > 0 
and 6Q^{t) > 0; when dC^{t) > 0, we have 5Q^{t) = 6Zn{t); and when dC^{t) > 0, we have 
Q?(t-) ~ z^lt-) ■ Hence, when 0 < Zn{t—) < Q^{t—), we have Zn{t) < Qn{t)- Moreover, the 
number of order arrivals for any given time horizon is finite with probability 1. □ 


This lemma, though simple, turns out to play an important role to ensure that fluid limits of 
order positions and related queues are well defined after rescaling. That is, we can extend the 
definition of and Z„ for any time t > 0. 
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For simplicity, for the rest of the paper we will use Q“, and Z„ instead of and 

Z„, defined on t > 0. The dynamics of the trnncated processes conld be described in the following 
matrix form. 


d 


Quit) 


/ 1 -1 -1 0 0 0 \ 

^ . 1-1-1 lQ5^(t_)>0,Qi(t-)>0,Z„(t-)>0 • d^n{t). (2.5) 

0 -1 -|||n 0 0 0 / 


The modified processes coincide with the original processes before hitting zero, which implies 
It<r„ = lQ-(t-)>0,Q^d-)>0,Z„(t-)>0- 

In order to establish the fluid limit for the joint process (Q(^, Q“ and Z„), we see that it is fairly 
standard to establish the limit process for from classical probability theory where various 

forms of functional strong law of large numbers exist. However, checking Eqn. (j2.5l) for Zn{t), we 
see that in order to pass from the fluid limit for to that for Z„(t), we effectively need to pass the 
convergence relation between some cadlag processes (X„, Y„) to (X, Y) in the Skorokhod topology 
to the convergence relation between f XndVn to f XdY. That is, consider a sequence of stochastic 
processes {X„}„>i defined by a sequence of SDEs 

Xn{t) = Unit) + f FniXn, S-)dYn(s), (2.6) 

where {U„}„>i, {Y„}„>i are two sequences of stochastic processes and {Fn}n>i is a sequence of 
functionals. Now, suppose that {U„, Y„, Fn}n>i converges to {U, Y, F} in some way. Then, would 
the sequence of the solutions to Eqn. (12.61) converge to the solntion to 

Xit) = U{t)+ f F{X,s-)dY{s)7 
Jo 

It turns out that such a convergence relation is delicate and can easily fail, as shown by the 
following simple example. 


Example 7 . Let {Xj}j>i be a seqnence of identically distribnted random variables taking valnes 
in {—1,1} snch that 

P(Xi = 1) = P(Xi = -1) = P(Xi+i = 1 I Xi = 1) = P(Xi+i = -1 I X, = -1) = ^ for i > 1. 

Define S'n(t) = Xi- Note that {Xi}j>i is a strictly stationary sequence, with mean zero 

and is a Markov Chain with finite state space { — 1,1}. Since each entry of the transition probability 
matrix is strictly between 0 and 1, the seqnence {Xi}j>i is -^-mixing, see e.g., [36]. Note that if:- 
mixing implies (/(-mixing, see e.g., m- By stationarity, 


lim —E 



= a 


CX) 

E[X2]+2^E[XiXi+i]. 

i=\ 


We can compnte by indnction that for any i > 1, 

E[XiW+i] = ^E[XiW] + iE[Xi(-Xi)] = ^E[XiW] = 







Therefore, we have 1 + ^ = 3. For strictly stationary centered (^-mixing sequence with 


E 


(ELi 


-)• oo as n —)• oo and E[|Xip+^] < oo for some 5 > 0, the invariance principle holds, 

see e.g., m, i-e., S„ converges to uB. Hence converges to \/3B. Now define a sequence of 
SDE’s dYn{t) = Yn{t)dSn{t) with Tn(0) = 1. Clearly, since Xi S {±1} and \Xi\ < 1, for sufficiently 
large n. 


\nt\ 






■v/n ^ — g y/n ‘ 


where |e„| < where C > 0 is a constant. Hence, converges to the limiting process described 
by exp{\/3.B(f) — |}, as n ^ oo. However, the solution to dY{t) = Y{t)d{'/3B{t)) with y(0) = 1 
is given by Y{t) = ex.p{y/3B{t) — 


Nevertheless, under proper conditions as specified in Assumptions [H [2l [5l one can establish the 
desired convergence relation. Such assumptions prove to be sufficient using a result of Kurtz and 
Protter @01 Theorem 5.4]. For sake of completeness, we present this result next, along with the 
technical conditions required for the convergence. 


2.3 Detour: Convergence of stochastic processes by Kurtz and Protter | I40] 

Define hs{r) : [0,oo) ^ [0, oo) by hs{r) = (1 — (5/r)+. Define Js : D^[0,oo) —)• £>”^[0,00) by 

^ hs{\x{s) - x(s-)|)(x(s) - x(s-)). 

S<t 


Let Yn be a sequence of stochastic processes adapted to T). Dehne Y^ = Yn — Js{Yn). Let 
Y^ = be a decomposition of Y^ into an T)-local martingale and a process with finite 

variation. 

Condition 8. For each a > 0, there exist stopping times r" such that P{r“ < 1} < 1/a and 
sup„ E[[M^]t< 7 -a + T{Af^)t<T^] < oo, where denotes the total quadratic variation of 

up to time r", and T{Aff)t<T^ denotes the total variation of Af^ up to time r“. 

Let Ti [0, oo) denote the collection of non-decreasing mappings A of [0, oo) to [0, oo) (in particular, 
A(0) = 0) such that \{h + t) — X{t) < h for all t,h > 0. Let be the space of real-valued k x m 
matrices, and oo) be the space of cadlag functions from [0, oo) to Assume that there 

exist mappings Gn,G : D^[0,oo) x ri[0, oo) — )■ D^km[0,oo) such that F/ o A = G„(x o A, A) and 
F{x) o A = G{x o A, A) for (x, A) G D^[0, oo) x ri[0, oo). 

Condition 9. (i) For each compact subset T-L C D^[0, oo) and t > 0, sup^^, supg<i \Gn{x, A, s) — 
G{x, A, s)| ^ 0; 

(a) For {{xn, A”)} G D^[0, oo) xTi [0, oo), sup 5 <^ |x„(s)—x(s)| —)• 0 and sup 5 <^ |A”(s)—A(s)| —)• 0 
for each t > 0 implies supg<j \G{xn-, A”, s) — G{x, A, s)| —)■ 0. 

Theorem 10. Suppose that (Un,X„,Y„) satisfies 

Xn{t) = Unit) + f Fn{Xn,S-)dYnis), 

Jo 
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(U„,Y„) => (U,Y) in the Skorokhod topology, and that {Y„} satisfies Condition for some 
0 < 6 < oo. Assume that {Fn} and F have representations in terms of {Gn} and G satisfying 
Condition\^ If there exists a global solution X of 

dX{t) = U{t)+ f F{X,s-)dY{s), 

Jo 

and the local uniqueness holds, then 

(U„,X„,YO^(U,X,Y). 

2.4 Fluid limit for order positions and related qneues 

We are now ready to establish our first result. 

Theorem 11. Given Assumptions\^\^\^ and\^ suppose there exist constants , q°‘, and z such 
that 


{Qii0),QU0),Zn{0))^iq\q^,z). 

Then, for any T > 0, Eqn. (j2.5p 

(Q^, Q“, Z„) ^ (Q^ Q“, Z) in {D^[0, T],J,), 
where (Q',Q“,z) is given by 

Q^{t) = q^ - A r), 

Q“(t) = - Au“(t Ar), 


and for t < t. 


dZ{t) 

dt 


Here r = min{T“, t°, r^} with 


= -A W + Y 


Au“ ’ 


r.3 Z{t-) 




Z{0) = z. 


r = 


Xvf>^ 


and 


( 11 ^ 


6(1 — e “6) 


c^{-l,0}, 

c = -1, 
c = 0. 

Moreover, if > 0, > 0, and > q^ jv^, then rf —>■ a.s. as n ^ oo, where 

a=\V^, b = q'-l(\V\ c=-|5. 


Mog(- + l 


= -V^ + + V^, = -V^ + V^ + V^. 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


( 2 . 12 ) 

(2.13) 
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Proof. Note that Eqns. (|2.7p . (12.8p . (12.91) satisfy the following SDE’s 



(Q\t)\ 

( 

d 

Q^t) 

= 


\z{t) ) 

\ 


/ 1 -1 -1 0 0 0 


0 0 


0 




0 -1 0 0 0 


^ ^ ^ ^Q‘^it-)>o,Q>>it-)>o,z{t-)>o^^dt; (2.14) 


(Q'’(0),Q“(0),Z(0))=(g^g^^). 

Therefore, we first show the convergence to Eqn. (j2.14p . Now, set = '^n, Xn = (Qn; Qn; 
and 


Fn{x,s-) = F{x,s-) = 


/I -1 
0 0 


-10 00 \ 
0 1 - 1-1 


0 - 1 


a;^(s—) 


0 


0 0 


)> 0 - 


/ 


In order to apply Theorem 1101 we need to decompose Yn- Now take 5 = oo, define the filtrations 

xr ■= o-({^(s)}o<s<nt,{f^i}i<j<Ar(no) ^nd Qi := cr{{^k}i<k<i), 


1 ^ 
Mn{t) = - V I 

n 

i=l 


and 


Ar,{t)=Yn{t)-M^{t). 

We will show that M„ is a martingale with respect to Ff and {l^}n>i satisfies Condition [8] in 
Theorem [TOl 

For s G [0, t), it is easy to see that Assnmption [3] implies that E[T^j | Ff i_i ] = E[T^* | Qi-i]. 

n Yk=l 


Thns 


E 


E 


= E 
= E 


I gi_i] \F^ n (iV(ns) < i) 

L n 2^k = l 

V^\F^n{N{ns) <i) 


F^ n {N{ns) < i) 


Meanwhile, Ff i n {N{ns) > i) F F^ n {N{ns) > i). Thus 

n Yk = l Fk 


E 


E 


1^ 


i I Qi—1 


= E 
= E 
= E 


E 




F^ n {N{ns) > i) 

F^ n (iV(ns) > i) 

I Qi-i n (iV(ns) > i) 


"^^\Fl. 
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Moreover, E 
Therefore, 


i\T^ n {N{ns) > i)l = since ^j is measurable with respect to (A^(ns) > i). 


E [Mn{t)\r;] = E 


'N{nt) 

V ^ 

i - E[T^ i Gi-i] 

-j^n 

i=l 

n 

'^S 


N{ns) 


^ ^ (^E I n (iV(ns) > i) 


-E 


i=l 


+ -E 
n 

N{ns) 


N{nt) 
i=N {ns) + l 


[e[^, I n (iV(ns) > i)] 
n (A^(ns) < i) 




I ^*-1 ^ ^ *)]) 

i=l 

+ ^Xn{t - s)(^E[t^i|j-," n (iV( ns) < i)] — E [e[ 1 ;^, I g,_i]|j-;n(iv( 77.5) < i) 

N{ns) 

-Y I ^*-1 ^ ^ *)]) = 


And E|M„(t)| < 00 follows directly from Assumption[2j Hence it follows that M„(t) is a martingale. 
The quadratic variance of Mn{t) is as follows; 


E[[v»i<] = 


J=1 


A - E 




^AE 

i=i 
6 


Vj] -2 VJe 




+ E 


n ^ 

j=i 


2 


< 


hift-i 


j=i 


since 


E 


V/E 




= E 


E 


V/E 




Qi-1 


= E E 


v/lQi-, 


Thus E[[M„] 4 ] is bounded uniformly in n since is square-integrable. Let [T{An)]t denote the 
total variation of An up to time t. Then E [[T(A„)]j] is also uniformly bounded in n, as 


E[[r(AO]i] =t^AE E[vi\gi.i 
i=i 
6 

= t^AE|H/| < 00 , 
i=i 


6 

< t J]AE[E[|t^^||6;i_i 
i=i 


(2.15) 

(2.16) 
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where the inequality in Eqn. (|2.15l) uses the Jensen’s inequality for conditional expectations and 
Eqn. (I2.16h follows from the square-integrability assumption. Thus, Yn satishes Condition [ 8 ] with 
r" = a + 1. Moreover, taking Gnix o e,e) = Fn{x) = F{x), it is easy to see that Condition [9] is 
satished according to [40]. 

Now = 0 when t as given in Eqn. (I2.10p : > 0 if — E® < 0; otherwise 

Q°‘{t) never hits zero in which case dehne r“ = oo. The case for is similar. 

It remains to find the unique solution for the limit Eqn. ()2.14l) . The equation for Z[t) when 
Z{t—) > 0 is a hrst-order linear ODE with the solution 


m 


- + At)) + 

(a log (6 — {t A r)) + z/b 
{z + ab)e~*^^ — ab 


\ ^ I + cj \b + c{t At) 
— alogb) ■ {b — (t A r)) 


l/c 




c = - 1 , 


(2.17) 


c = 0 . 


From the solution, we can solve explicitly as given in Eqn. ()2.1ip . Note that the expression for 
Z{t) may not be monotonic and there might be multiple roots when c 7 ^ 0. Nevertheless, it is easy 
to check that the solution given in Eqn. (|2.11l) is the smallest positive root. For instance, when 

c ^ {— 1 , 0 }, there are two roots —b/c and + b^ ^ when c = — 1 , 

there are two roots b and 6(1 — e~^). More computations confirm that indeed the smallest positive 
roots are + 6 ^ 5 i/(c+i)g-i _ fQj. ^ ^ {— 1 , 0 } and Y = 6(1 — e“^) for c = — 1 . 

Moreover, from these calculations. Therefore r = min{r“, r^j is well defined and hnite. □ 

The following hgure illustrates the fluid limits of (Q^(t), (5“(t), Z(t)) with (5^(0) = (5“(0) = 
Z(0) = 100, X = 1,V^ = V^ = l,V‘^ = 0.6, = 0 . 8 , = 0.7, E® = 0.8. 


Q^,Q^,Z against time 



Figure 2: Illustration of the fluid hmit(Q^(f), Q“(f), Z(f)). 
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3 Fluctuation analysis 


The fluid limits in the previous section are essentially functional strong law of large numbers, and 
may well be regarded as the “first order” approximation for order positions and related queues. 
In this section, we will proceed to obtain a “second order” approximation for these processes. We 
will first derive appropriate diffusion limits for the queues, and then analyze how these processes 
“fluctuate” around their corresponding fluid limits. In addition, we will also apply the large de¬ 
viation principles to compute the probability of the rare events that these processes deviate from 
their fluid limits. 


3.1 Diffusion limits for the best bid and best ask queues 

We will adopt the same notation for the order arrival processes as in the previous section. However, 
we will need stronger assumptions for the diffusion limit analysis. 

There is rich literature on multivariate Central Limit Theorems (CLTs) under some mixing 
conditions, e.g.. Tone [48]. However, these are not functional CLTs (FCLTs) with mixing conditions. 
In the literature of limit theorems for associated random fields, FCLTs are derived under some weak 
dependence conditions with explicit formulas for asymptotic covariance of the limit process. Here, 
to establish FCLTs for {T^i} i>i, we will follow as in m- Readers can find more details in the 
framework of Bulinski and Shashkin m Chapter 5, Theorem 1.5]. 

Assumption 12. {N{i,i -|- l]}i£z is a stationary and ergodic sequence, with A ;= E[A(0, Ij] < oo, 
and 

OO 

^||E[iV(0,l]-A| JT-]|| 2 <oo, (3.1) 

n=\ 

where ||T ||2 = (E[y^])^/^ and '■= -|- l],i < —n). 

Assumption 13. Let n G N and A4(n) denote the elass of real-valued bounded coordinate-wise 
non-deereasing Borel funetions on M"'. Let |/| denote the cardinality of I when I is a set, and || • || 
denote the L°°-norm. Let {^i}i>i be a stationary sequence ofM.^-valued random vectors and for 
any finite set / C N, J C N, and any f,gG A4(6|/|), one has 

Covif(^i),g(^j))>0. 

Moreover, for 1 < j < 6, 

OO 

= Var{vi) + 2 ^ Cov{Vi,V^) < oo. 
i=2 


Remark. Note that an i.i.d. sequence {Vi}i>i clearly satisfies the above assumption if 1/ i 
is square-integrable. It is not difficult to see that Assumption [12] implies Assumption [U and 
Assumption [13] implies Assumption [2] In particular, Theorem 1111 holds under Assumptions [12] and 

m ^ _ 

With these assumptions, we can define the centered and scaled net order flow = ('I'„(t) ,t> 
0 ) by 



' N{nt) 






2=1 




(3.2) 
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Here 


^ = (V\ 1 < j < 6 ) = (E[h/], 1 < i < 6), 
is the mean vector of order sizes. 

Next, define and R“, the time rescaled queue length for the best bid and best ask respec¬ 
tively, by 

dRiit) = + XVH) - di^lit) + XVH) - di^lit) + XVH), 

dRi{t) = d(<(t) + XVH) - di'I’lit) + XVh) - d(<(t) + XV%. 


The definition of the above equations is intuitive just as their fluid limit counterparts. The only 
modification here is that the drift terms is added back to the dynamics of the queue lengths because 
^ has been re-centered. The equations can also be written in a more compact matrix form. 



with the linear transformation matrix 

/I -1 -1 0 0 0\ 

“Vo 0 0 1 -1 -1 y ■ 


(3.3) 


(3.4) 


However, Eqn. ()3.3p may not be well defined, unless R^it) > 0 and R“(t) > 0. As in the fluid limit 
analysis, one may truncate the process at the time when one of the queues vanishes. That is, define 


C = : Kit) < 0}, = inf{t : R^(t) < 0}, Ln = mf{C, (3-5) 

and define the truncated process (R^,R“) by 



Now, we will show 

Theorem 14. Given Assumptions\^ liM and MA for any T > 0, 


• We have 

^ o Ae - HurfWi o Ae in {D^[0,T], R). (3.7) 

Here Wi is a standard scalar Brownian motion, is given by Eqn. ()3.1ip . ^ is a standard 
six-dimensional Brownian motion independent o/Wi, o denotes the composition of functions, 
and Ti is given by TjTj'^ = [ajk) with 


- / 

for 

j = k, 

^jk — \ 

y Pj^k'^j^k 

for 

j^k, 


(3.8) 
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and 


v] = VaHyi) + 2 ^ Cov{vi, 




i=2 


/ ^ 

p.^ = ^{coviyl,v,^) + Y. 


(3.9) 


'l,Vn + Cov{VtV^ 


i=2 


That is, < J < 6 ) zs a six-dimensional Brownian motion with zero drift and 

varianee-eovarianee matrix (AS^S + Xv^V ■ V'^). 


V (^n(0),-Rn(0)) ^ (g^g“), then for any T > 0, 


R“ 


R® 

R“ 


in {D^[0,T],Ji). 


Here, the diffusion limit process (R^, R“)^ up to the first hitting time of the boundary is a two- 
dimensional Brownian motion with drift 'jt and the varianee-eovarianee matrix as 


ft := (/xi,//2)'^ = • y 

Proof. First, define N„ by 


and 


aa'^ := A • (AS^S + Xv^V • V^} • A^. 


(3.10) 


Nn(t) = 


N{nt) — nXt 


n 


Now recall the FCLT from [51 Page 197]. For a stationary, ergodic, and mean-zero sequence 
{Xn)nez, that satisfies ||E[Xo | XZff ]\\2 < oo, we have Wi{-) on (iA[ 0 ,r], Ji) 

with = E[Xq] -|- 2^^^E[XoX„] < oo, where Wi is a standard one-dimensional Brownian 
motion. Since the sequence {N{i,i -\- Ijjigz satisfies Assumption [T^ 


fV^n-J - A[n-J 


n 


VdWif 


in (iA[0, T], Ji) as n —oo, where 


vl = E[(iV(0,1] - A)2] + 2 J]E[(iV(0,1] - A)(iV(j, j + 1] - A)] < oo. 

1=1 

Next, for any e > 0 and n sufficiently large, 

X^lns] - X [nsj Nns - Ares 


P sup 
V 0<s<T 


> e 


< 


ytn y/fi 

N[k, A: -|- 1] > e^/n — x}j 

< ([reTj + 1)P(1V[0,1] > ev^ - A) 

[reTj -M f 112 


max 
0<k< [nTJ ,/cGZ 


< 


/ 

Jn 


- A)^ 7Ar[0,l]>ev^-A 


iV[0, Ij^dP ^ 0, 


(3.11) 
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as n —>■ oo. Hence, => I’dWi in {D[0,T], Ji) as n ^ oo. 

Moreover, thanks to |16[ Theorem 2], Assumption 1131 implies 

in (Z)6[0,r],Ji), 

where ^ is a standard six-dimensional Brownian motion and S is a 6 x 6 matrix representing the 
covariance scale of the limit process. Fnrthermore, the expression of S by Eqn. (13.8p and Eqn. (13.9p 
can be explicitly compnted following m Theorem 2]. 

Now, by Assumption [3l the joint convergence is guaranteed by [53l Theorem 11.4.4], i.e., 

in (^^[O, T], Ji). 

Moreover, by m Corollary 13.3.2], we see 

EudWioAe in (T»®[0, T], Ji). 


To establish the second part of the theorem, it is clear that the limiting process would satisfy 


d 




(3.12) 


with 

= inf{t : < 0}, = inf{t : < 0}, l = min{t“, (3.13) 

We now show that 

(R^,R“)^(R',R“) in (Zl2[o,r],Ji). (3.14) 

According to the Cramer-Wold device, it is equivalent to showing that for any {a,/3) G 


aR^ +/3R“ ^ oR'’+/3R“ in {D‘^[0,T], R). (3.15) 

Since ^ ^ in (Z1^[0, T], Ji), by the Cramer-Wold device again, 

in {D^[0,T], R). 


By dehnition, it is easy to see that 

aR^(t) -(- f3R^{t) = (a, /3) ■ A ■ n{t A in) + V {t A in)^ + aq^ + ldq°‘■ 

Since the truncation function is continuous, by the continuous-mapping theorem, it asserts that 
Eqn. (I3.15P holds and the desired convergence follows. 

Moreover, because Ee is deterministic and is a constant, we have the convergence in 

Eqn. (I3.15p . as well as the convergence in Eqn. (|3.14l) . Note that in-,'n>l and t are hrst passage 
times, by [53l Theorem 13.6.5], 

{in, Ri{in-),Rl{in-)) ^ (i, R\i-), R\i-))- □ 
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Remark. Theorem 1141 holds without Assumption [3l as long as is guaranteed to con¬ 

verge jointly. 


3.2 Fluctuation analysis of queues and order positions 

Based on the diffnsion and fluid limit analysis for the order position and related queues, one may 
consider fluctuations of order positions and related queues around their perspective fluid limits. 

Theorem 15. Given Assum,ptions^ f5l liHL and MA we have 

/ Qn - \ / Trl - _ ^3 \ 

Vni Q“-Q“ ^ 1 in {D^[0 ,t),Ji) 

VZn-Z/ V Y J 


as n —>■ oo. Here (QtQ“,ZO, (Q^Q^ Z) are given in Eqn. (12.21) and Theorem Mli , 1 < j < 6) 
is given in Eqn. (1321), and Y satisfies 


dY(t) = 


Z(t)(Ti(t)-T2(t)-T3(t)) 



(3.16) 


with Y (0) = 0. 

Proof. Given Assumptions El El [121 and fT31 we have from Theorem flTl 




' N{n-) 

E 

2=1 


- XnVe 




in (T»®[0 ,r), Ji). 


Hence, we have the following convergence in (D[0,t), Ji), 

- Q'’) ^ _ ^3^ 

Vn{Qn - Q“) ^ - T'®. 


Since Theorem 1111 holds under Assumptions [T2] and 1131 we now use the dynamics of Zn{t) in 
Eqn. (12.5p and Z{t) in Theorem 1111 and get 


d{Zn{t)-Z{t)) = -diCf{t)-C^{t))- 


2t.\ r^2/.\\ ^n{t ) ,^3u\ I 


QUt-) 


dCffit) + 


QKt-) 


dC\t) 


= -d{Cl{t) - C\t)) - ^E^MClit) - C\t)) + 


We can rewrite this as 

d{Zn{t) - Z{t)) + 

Zn{s-) 


QUt- 

Zn{t-) - Z{t-) 
Q\t-) 


Z{t-) Zn(t-) 

LQHi-) Q’kit-) 


dC^it). 


dC\t) = dXn{t), 


X^{t) = -{Cl{t)-C\t))- [ 

Jo 


0 Qnis~) 


dwUs) - c\s)) + f 

Jo 


Q\s-)Q\{s-) 
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Now, 


^ X „ ^ -^2 _ 




'{s)+ [ 

Jo 


lo Q'^is-) ^' Jo ms-w 

As the limit processes ^ and Q^, Q“ are continuous, this could be changed into 


W^ds 


^ -^2 _ / ^ 

Jo h/ i'Sj Jo 


Zis) , r Zis){^Hs) - ^Hs) -IM 




Hence, 

where Y satisfies Eqn. (j3.16l) . 


v^(Z„-Z) ^ Y, 


□ 


3.3 Large deviations 

In addition to the fluctuation analysis in the previous section, one can further study the proba¬ 
bility of the rare events that the scaled process , Qni'^)) deviates away from its fluid limit. 

Informally, we are interested in the probability P(((5^(f), Q“ (f)) ~ f°'{t)),0 < t < T) as 

n —)• oo, where (/^(t),/“(t)) is a given pair of functions that can be different from the fluid limit 
{Q^{t),Q-{t)). 

Recall that a sequence (Pn)neN of probability measures on a topological space X satisfies the 
large deviation principle with rate function X : X —S' M if X is non-negative, lower semi-continuous 
and for any measurable set A, we have 

— inf X(x) < liminf — logXn(^) < limsup — logP„(A) < — inf I{x). 
xeA° n n^oo n xeA 

The rate function is said to be good if the level set {x \ I{x) < a} is compact for any a > 0. Here, 
A" is the interior of A and A is its closure. Finally, the contraction principle in large deviation says 
that if Pn satisfies a large deviation principle on X with rate function X(x) and X : A —)■ Y is a 
continuous map, then the probability measures Qn ■= PnF~^ satishes a large deviation principle 
on Y with rate function I{y) = \iiix\p(^^'^^yZ{x). Interested readers are referred to the standard 
references by Dembo and Zeitouni |24] and Varadhan [49] for the general theory of large deviations 
and its applications. 

Recall that under Assumptions [Il[2| and (H we had a FLLN result for ((5^(f),<5“(t)) and under 
Assumptions US [HI andO we had a FCLT result for {Q^{t), Q'^{t)). It is natural to replace As¬ 
sumptions [H [2] by some stronger assumptions to obtain a large deviations result for (Q^(t), Q“(f)). 
We will see that by assuming that (I^i)i6N and {N{i) — N{i — I))j6N satisfy the following Assump¬ 
tions dsi and [13 in addition to Assumption [3l by a large deviation result of Bryc and Dembo |14j , 
we will have the large deviations for , Q'}^{t)). 

Assumption 16. Let be a sequence of stationary -valued random vectors with the a- 

algebra defined as a{Xi,m < i < £)■ For every C < oo, there is a nondecreasing sequence 
£{n) € N with n{n+i) ^ 

sup {P(A)P(H) - nB)\AGF^\BG G n} < 

sup |P(A n H) - e^(’")p(A)P(H) \AeFlf\Be G n} < 
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Assumption 1161 holds under the hypermixing condition in [241 Section 6.4], under the ^-mixing 
condition of Bryc |13[ (1.10),(1.12)], and nnder the hyperexponential a-mixing rate for stationary 
processes of Bryc and Dembo [m Proposition 2]. Therefore, if (l^j)ieN and {N{i) - N{i - l))i 6 N 
satisfy Assumption [TBl then Assumptions [T] and [2] are satisfied. It is also clear that Assumption fTHl 
holds if Aj’s are m-dependent. 

In order to have the large deviations result, we also need to assume that (1^ i)i£N and {N{i) — 
N{i — l))j 6 N satisfy the following condition: 


Assumption 17. For all 0 < ’y,R < oo, 


9r{i) ■= 


sup — log E 

fc,mSN,fcS[0,-Rm] ^ 


e7llE 


fc+m 
i=k + l 


Xi 


< oo, 


and A := sup..^ limsup^^oo ^ ^9r{i) < oo. 

Note Assumption [T7] is trivially satisfied if Xj’s are bounded. If Xj’s are i.i.d. random variables, 
Assnmption [T7] rednces to the finiteness of the moment generating function of Xi, which is a 
standard assnmption for Mogulskii’s theorem l |24l Theorem 5.1.2]). Therefore, Assnmption 1171 is 
a natural assumption for large deviations. 

Under Assumption [16] and Assnmption (TT] Dembo and Zajic |23j proved a sample path large 
deviation principle for P(^ '^i=i ^ ') (For ease of reference, we list it in Appendixj^as Theorem 

[38|) . From this, we can show the following. 

Lemma 18. Let both (T^i)i 6 N {N{i) — N{i — l))igH satisfy Assumvtion[T6\ and Assumption[T7\ 
and let Assumption\^hold. Then, for any T > 0, P(C'„(t) G •) satisfies a large deviation principle 
on Loo[0,T] with the good rate function 


I{f) = inf [Iv{g) + lN(,h)], 

heACf [ 0 ,r] ,g£ACo [0,oo) 
g{hlt))=f{t),0<t<T 


with the convention that inf 0 = oo and 

poo 

Iv{9) = / ^v{ 9 '{x))dx, 

Jo 

if g € oo) and Iv{9) = oo otherwise, where 


Ay(x) := sup {0 • X — r\/(0)} , ry(0) := lim — logE 

n—)-oo Tl 




(3.17) 


(3.18) 


and 


if h £ AC(j’[0,T] and Inih) 


Ixih) = [ AAr(/i'(x))dx, 

Jo 

oo otherwise, where 


An{x) 


sup {0 • X - rAr(6»)} , rAr(6') 

0eR6 


lim — log E 

n—>cx) Tl 



(3.19) 
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Proof. Under Assumption! 161 and Assumption [T71 by Theorem 138! in Appendix|3 P(^ S[=i Vi ^ ') 
satisfies a large deviation principle on Loo[0, M] with the good rate function 

j-M 

Ivif) = / Av{fix))dx, 

Jo 

if / G ACq[0,M] and Iv{f) = oo otherwise, where Ay(x) and ry(0) are given by Eqn. (I3.18|) and 
P(iA^(n-) G •) satisfies a large deviation principle on Loo[ 0 ,T] with the good rate function 


InU) = [ ANif{x))dx, 
Jo 


if / G ACq[Q, T] and InH) = cc otherwise, where An{x) and PNid) are given by Eqn. (I3.19p . Since 
and Nt are independent, P(^ I]i=i i G •, ^N{n-) G •) satisfies a large deviation principle 
on Loo[0,M] X Loo[0,T] with the good rate function Iv{-) + 

We claim that the following superexponential estimate holds: 


limsuplimsup — logP (A'(n) > nM) = —oo. 
M—>-oo n^oo n 


(3.20) 


Indeed, for any 7 > 0, by Chebychev’s inequality. 


P(iV(n) > nM) < 




Therefore, 


lim sup — log P {N (n) > nM) < 

n—)-oo ^ 


—7 + lim sup — log E 

n—)-oo 


^'yN{n) 


(3.21) 


From Assumption [T71 sup..y>Q limsup^^oo ^ logE < 00 . Hence, by letting 7 —)• 00 in 

Eqn. (j3.2ip . we have Eqn. ()3.2np . 

For any closed set C G Loo[0, T], 


1 

lim sup — log P — Vi £ C 

„ > r) * 77, ‘ ^ 


n^oo ^ 


2=1 


iV(n-) 


= lim sup lim sup — log P ( — i £ C, —N inT) < 

n T 71 \ 71 71 


2 = 1 


M 


[ivig) + iN{h)] 


M^oo n—>-oo n 

= — inf inf 

MeN f&C 

h&ACf{0,T] ,g&ACo [0,M] 
9(h{t))=f{t)A<t<T 
h{T)<M 

= -inf inf [Iv{g) + lN{h)], 

fee heAC+ [0,T] ,geACo [0,oo) 
g(h(t))=f{t)fl<t<T 


(3.22) 

(3.23) 


where Eqn. (I3.22p follows from Eqn. ^3.201) and Eqn. p3.23p follows from the contraction principle. 
The contraction principle applies here since for h{t) = ^N{nt) and g{t) = ^Yli=i i have 
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n ^~ 9 {h{t)) and moreover, the map {g, h) i—)• goh is continuous since for any two functions 

Fn,Gn —>■ -F, G in uniform topology and that are absolutely continuous, we have sup^ |F„(Gn(t)) — 
F{G{t))\ < supt \Fn{Gn{t)) - F{Gn(t))\ + sup^ \F{Gn{t)) - F{G{t))\ 0 as n -)► oo. 

For any open set G G Loo[0,T], 


lim inf — log P 

n^oo Tl 



N(n-) 


> lim inf — log P ( — € G, —N{nT) < 

n^cxD n \ n n 


M 


2 = 1 


inf 

f&G 

heAC^ [o,r] ,geACo [o,m] 

g{h{t))=f{t),0<t<T 

h{T)<M 


[ivig) + lN{h)]. 


Since it holds for any M G N, the lower bound is proved. □ 

Moreover, by the contraction principle. 

Theorem 19. Under the same assumptions as in Lemma\JR P((Q^(t), Q® (t)) G •) satisfies a large 
deviation principle on L°°[0, oo) with the rate function 

/(/',/“)= inf m, 

(p^yf 

where !{■) is defined in Lemma \JR Qf is the set consists of absolutely continuous functions 4>{t) 
starting at 0 that satisfy 


d{fHt),r{t)f 


1 - 1-10 0 0 
0 0 0 1 -1 -1 


d(j)(t), 


with the initial condition (/^(O),/“(O)) = {q^,q°‘)- Otherwise I{f) = oo. 


Proof. Since P(Cn(t) G •) satisfies a large deviation principle on L°°[0,oo) with the rate function 
Fif)), it follows that F{{Q^{t),Qf^{t)) G •) satisfies a large deviation principle on L°°[0, oo) with the 
rate function 


/(/) 


inf F((j)) 


where Qf is the set of absolutely continuous functions 4>{t) = 
satisfy 


(^ / 1 -1 -1 0 0 
VO 0 0 1-1 


{fd (t), 1 < j < Q) starting at 0 that 

^ #(t), 


with the initial condition (/^(O),/“(O)) = {q^,q°‘). It is clear that 


/(t) = q’' + (j)^{t) - 4>‘^{t) - #(t), 

r{t)=q‘^ + ^\t)-^Ht)-^Ht), 
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and the mapping 4> {f^, /“) is continuous, since it is easy to check that if 

4>n{t) ■■= (p{t) = 

in the L°° norm, then (/n(t),/n(i)) “S' (/^(^))/“(^)) in the L°° norm. Since the mapping (j) i-)- 
(/^,/“) is continuous, the large deviation principle follows from the contraction principle. □ 

Let us now consider a special case of Theorem [19) 


Corollary 20. Assume that N{t) is a standard Poisson process with intensity X independent of 
the i.i.d. random vectors in such that < oo for any 0 E M®. Then, the rate function 

I{f) in Eqn. (I3.17P in Lemma {TR has an alternative expression 


lif) = / A{nt))dt, 


(3.24) 


for any f E ^Co[0,oo), the spaee of absolutely continuous functions starting at 0 and = +oo 
otherwise, where 

A{x) := sup io-x — A(E[e®'^i] — 1)| • 


Proof. First, notice that when Nt is a standard Poisson process with intensity A, independent of 
i.i.d. random vectors then, N{i) — N{i — 1) is a sequence of i.i.d. Poisson random variables 
with parameter A and therefore both and {N{i) — N{i — l))igN satisfy Assumption [Tbl and 

Assumption E] is also satisfied. Under the assumption, E[e®'^i] < oo for any 0 E M® and 
jgjg6»(Af(i)-Af(i-i))j _ gA(e®-i) ^ 0 E M. Therefore, both ( 1 ^ 1 ) 161 !} and {N{i) — N{i — 1)) 

satisfy Assumption (TT) 

By Lemma dS) 


moreover, 
ieN 


pi poo 

Iv{g) + lN{h)= / Avig'{t))dt+ / AN{h'{t))dt, 
Jo Jo 


where 


and 


Ay(x) = sup • X — logE | 

0 eR® L J J 


An{x) = xlog - X + A. 
Since fft) = g{h{t)), we have f'{t) = g'{h{f))h'{t) and 


J Av{g'{t))dt = j Av{g'{h{t))h'{t)dt = j Ay h'{t)dt. 


Therefore, 


inf {Iv{g) + lN{h)) 

heACf[0,T],geACo[0,cxy) 
g{h{t))=f{t),o<t<T 


= inf r (av 
h&ACf[0,T\Jo V 



h'{t) + h'{t) log 



h'{t) + X 


dt. 
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Now, 


= inf sup io ■ X — y log]E[e®'^M + y log (^') 
y e ^ VA/ 

= sup inf io ■ X — y loglEle^'^M + y log 
e y ^ VA 

= sup • X — A(E[e^''^^] “ 1)| • 


-y + A} 

-2/ + a} 


Therefore, Eqn. ()3.17l) reduces to Eqn. ()3.24p . 


□ 


4 Applications to LOB 

4.1 Examples 

Having established the fluid limit and the fluctuations of the queue lengths and order positions, 
we will give some examples of the order arrival process N{t) that satisfy the assumptions in our 
analysis. 

Example 21. [Poisson process] Let N{t) be a Poisson process with intensity A. Clearly assump¬ 
tions [U and m are satisfied. 


Example 22. 

with intensity 


[Hawkes process] Let N{t) be a Hawkes process [T2], i.e., a simple point process 


\{t) := A 



s)N{ds) , 


(4.1) 


at time t, where we assume that A : M>o —)• M"*" is an increasing function, a-Lipschitz, where 
h : M>o —)• M"*" is a decreasing function and h{t)tdt < oo. Under these assump¬ 
tions, there exists a stationary and ergodic Hawkes process satisfying the dynamics Eqn. m (see 
e.g., Bremaud and Massoulie m)- By the Ergodic theorem. 


A:=E[iV(0,l]], 

a.s. as t —>■ oo. Therefore, Assumption [T] is satisfied. It was proved in Zhu [55], that {A^(i, i-t-l]}igz 
satisfies Assumption [12] and hence '^dWii')-, on (H[0,T], Ji) as n —>■ oo. 

In the special case A( 2 :) = v + z, Eqn. ( 14 . 11 ) becomes 


A(t) — V -\- 



s)N{ds), 


which is the original self-exciting point process proposed by Hawkes m, where u > 0 and II^IIli < 
1. In this case. 


A 


V 



V 

{i-\\h\\L^r' 
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Example 23. [Cox process with shot noise intensity] Let N{t) be a Cox process with shot noise 
intensity (see for example 0). That is, N{t) is a simple point process with intensity at time t given 
by 


/ 


\{t) = u + / g{t — s)N{ds), 


where is a Poisson process with intensity p, g : M>o —>• K'*' is decreasing, jj^llii < oo, and 
tg{t)dt < oo. N{t) is stationary and ergodic and 


m 

t 


X := 


^ + pWoWl^ a.s.. 


as t —)■ oo. Therefore, Assumption [T] is satisfied. Moreover one can check that condition Eqn. (13.ip 
in Assumption [12] is satisfied. Indeed, by stationarity. 


||E[iV(0,1] - A| = ||IE[iV(n - l,n] - A| 


We have 


where 


therefore. 


E[A^(n - 1, n] - A I To"°°] = IE 


X{t)dt — A 


/ n—1 


T 7 —oo 
*^0 


\{t) = u + f g{t — s)N{ds)+ f g{t — s)N{ds), 

J —OO J 0 


rn pO pn pt 

E[A^(n — 1, n] — A I = j j g{t — s)N{ds)dt + p j / g{t — s)dsdt — p\\g\\i^i. 

Jn-i Jo 


By Minkowski’s inequality, 
||E[iV(n-l,n]-A| J-o-°°]|| 2 < 
Note that 


' n—1 J —oo 


pn pD pn pt 

/ / g{t - s)N{ds)dt +p / g{t - s)dsdt - pUgW^i 

Jn—lJ—co 2 Jn—lJo 

pn pt pn poo 

P / ait-s)dsdt-p\\g\\Li =p / gis)dsdt, 

Jn—lJo 2 Jn—lJt 


therefore. 


pn pt poo poo poo 

y] P / / ait - s)dsdt - p\\g\\Li = / gis)dsdt= / tgit)dt. 

n—i Jn—1 Jo 2 Jo Jt Jo 
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Furthermore, 


E 

72=1 


rn /*0 


' 71—1 J —OO 


g{t — s)N{ds)dt 


< 


g{n — 1 — s)N{ds) 


2 n=l 

J —OO 

OO 


= E\ 

n=l > 

/ / 9^ 

J —OO 

OO 

If 

J —OO 

OO 

^ Vp'^Vain 

71=1 

VI 

OO 

'^9{n 

71=1 


OO „o 


(n — 1 — s)ds 


(n — 1 — s)ds 


— 1 — s)ds + p tg{t)dt 

Jo 


CXD rtQ 

^ff(n - 1) + ^ / g{n-l- s)ds 

n=l n=l 


+ P 


roo 

/ tg{t)dt 
Jo 


< 


Vp 


5 ( 0 ) + Ibllii + / tg{t)dt 
Jo 


+ p tg{t)dt < OO. 

Jo 


Hence Assumption [12] is satisfied. 


Nri-Xn- 


—^^ VdWi{-) in (I7[0, T], Ji) as n ^ 00 , where 

+ plbllii + 


4.2 Probability of price increase and hitting times 

Given the diffusion limit to the queue lengths for the best bid and ask, we can also compute 
the distribution of the first hitting time i (defined in Eqn. ()3.13p l and the probability of price 
increase/decrease. Our results generalize those in [20] which correspond to the special case of zero 
drift. 

Given Theorem 1141 let us first parameterize a by 

/ p^ aip \ 

V 0 aa y’ 

and assume that —1 < p < 1. Next, denote ly the modified Bessel function of the first kind of 
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order v and t'n := revr/a, and define 


a := < 


TT + tan 


-1 




tan 


-1 




p > 0, 
p = 0, 
p<0, 


ro := 


{q’^/cTiy + (g“/fT2)2 - 2p{q’’/ai){q°-/a2) 


-1 / qr 


Oo '■= 

/ . 

(}'>/(Ti-pq“/cr 2 

Then according to Zhou [21], we have 


1 — 

-1 g°/o- 2 -\/l-p^ 


vr + tan ^ - 5-7 -— 

' q / <^2 


q’^/ai < pg“/cT2, 
q’^/ai = pg“/cJ 2 , 
q'^/ai > pq°‘/(J 2 . 


Corollary 24. Given Theorem I 4 and the initial state ( 9 ^,( 7 “), the distribution of the first hitting 
time i 


i>t) = —fiiq'’+i2q°-+ht 
at 


J^sin 


/ mrOi 


n=l 


\ a 


o\ -li r 

e 2 ‘ / sin 


pa 

/ sii 

Jo 


mrO 


a 


gn{0)d9, 


where 


g^(e) := 

ppia 2 - P 2 cri _ 12^2 ,7 I 7 

■, I 3 -h ptlt20'l<72 H--l-'iPi + '2P2, 


_ -pifT 2 + PP2<71 , _ 


ifaf , , loCTo 


(1 - p 2 )cjfcr 2 


(1 -P^) 0 ' 2 f^l 


U := ZlUi + pl2Cr2, k ■= /2<72\/1 - p^- 


Note that when ft > 0, it is possible to have = 00 ) > 0, meaning the measure P-^ might be 
a sub-probability measure, depending on the value of 7^. In this case, P-^(i. > t) actually includes 

P;^(i = 00). _ 

Moreover, based on the results in Iyengar [36] and Metzler 


Corollary 25. Given Theorem If and the initial state {q^,q°'), the probability of price decrease is 
given by 

rcxD roo 

= / / exp(K^(rcosa — -|-it“(rsina — 2 :“) — |’^|^t/ 2 )p(t,r)dr(it, 

Jo Jo 

'mT{a — 60 )' 


where 


g{t,r) = ^ e (^^+^o)/ 2 i nsin 

a^tr 

n=l 


a 


T 

^n-Kja [,17 


and it = [tJ^, k°'Y'= a ^{pi, ^ 2 )'^ and {z^, z°') = a That is, 


K 


<72P1 _ 

-fJippi -b - P‘^P2 


<72 r 


-cTipg^ -b CTIV^I - p‘^q°' 
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Similarly, when 'jt > 0 , with positive probability, we might have = oo and = oo. Therefore 
we compute here implicitly refers to < <■“, < oo) in that case. 

Note that both expressions for i and the probability of price decrease are semi-analytic. However, 
in the special case of 7 ^ = 0 , i.e., when Vi = ¥2 + H 3 and V 4 = + Vi^, they become analytic. 

Corollary 26. Given Theorem [771 and the initial state {q^,q^), if~ft = if, then 


P(r > t) 




n: odd 


1 

— sin 
n 


uttOq 

a 




Corollary 27. Given Theorem I4 and the initial state if~it = 0 , then the probability that 

the price decreases is 


Proof. 


P(i^ < i“) = 


/ 


^sin(7r0o/a) 


dr 


/ 


0 sin^(7r0o/a) + ((r/ro)'^/'^ + cos{tt9q/ a))"^ arg 
r°° sin(7r0o/a) d{r/rQy^°‘ 

0 sin^(7r0o/a) + + cos(7r0o/a))^ 

sin(7r6*o/a;) dx 6 q 

0 sin^(7r0o/a) + (a: + cos(7r0o/Q;))^ ^ ot 


□ 


4.3 Fluctuations of execution and hitting times 

In addition, we can study the fluctuations of the execution time rf. 
Proposition 28. Given Theorem 1 1 51 for any x (say, x <0), 


lim ¥{^/n{Tf — r^) > x) = P(y(r^) > ax) = 1 — ^ f — ^ 

n-^00 ycTyfr^j y 


where <I>(x) := 


= -!/T 2 


dy is the cumulative probability distribution function of a standard Gaus- 


-00 

sian random variable, with cry the variance of the Gaussian process Y{t). 


Proof. For any x < 0, 


nv^i<-T^)>x) 

- F (vi {Z. . (r- . ^)) > -Vi. (.» + ^)) . 


Note that 

lim \fnZ (t ^-\—= xZ'{t^), 

n^oo \ y/n) 

and for any t > 0, c 7 ^ — 1, Eqn. (I2.17P and Eqn. (12.111) lead to 


Z'(r^) 


ac 

1 + c 



^(1 + c)^ 



l + c 
c 


= —a. 
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Similarly, when c = —1, we have Z{t) = (alog(6 ~ ^) + f ~ a\ogb){h — t). Thus Z'{t) = —a — 
(olog(6 — i) + f — a log 6), and Z'{t^) = —a — ^alog(6e“'^) + f — a log 6^ = —a. Finally, recall 

that ^/n{Zn{t) — Z{t)) —)• Y{t) on {D[0,t^), Ji) as n —)• oo, hence the first equation. 

The second equation follows from Y(t) being a Gaussian process with zero mean and variance 
CTy, the latter of which can be computed explicitly, albeit in a messy form as in Appendix iBl □ 

Proposition 29. Given Theorem m with v^,v°‘ > 0 , for any x (say x <0), 

(^) 


lim — r“) > x) = 1 — 


q^Xv^ 


"011 + ■022 + 'f’SS - 2V’12 - 2'iAl3 + 2V’23 


X 


(4.2) 


where 4>(x) := /-oo ^ is the cumulative probability distribution function of a standard 

Gaussian random variable, and 


fjij := ^ SifcS.feA + V^V^vjX^ l<i,j <6. 


(4.3) 


k=l 


(li) 


lim P(v^(r“ - r“) > x) = 1 - $ 


q°‘Xv° 


■044 + V'SS + ^66 — 2-045 “ 2046 + 2056 


(4.4) 


Proof. Similar to the proof of the fluctuation of the execution time rf, we can show that, for any 
a: < 0 , 

lim P(v^« - r^) > x) = P((T'^ - ^2 _ ^ -(Q^)'(r^)x), 

n—^oo 

From the expression of in Eqns. (12.71) . (12.101) . and (13.7p . it is clear that = —q^ and 

the mean of (T'^ — — T'^)(t) is zero and the variance is 


(011 + 022 + 033 - 2012 - 2013 + 2023 )t- 


Therefore, 


q^Xv^ 1 

- X I 

011 + 022 + 033 - 2012 - 2013 + 2023 j ' 

Similarly, we can show that Eqn. (|4.4I) holds. □ 

Finally, we have the large deviations for the tails of the hitting time. Indeed, given Assumptions 
dl [2]and[3l we have the fluid limit in Theorem [TTl we see that —>• r := A r“. More generally, 
by replacing Assumptions [T] and [ 2 ] by the stronger Assumption fTHl we can use the large deviations 
result to study the tail probabilities of the hitting time as n goes to oo. Note that for any t > t, 

^{Tn >t)=F (^Qn{s) > 0, Q“(s) >0,0<s<t) = P (^Qn{s) > 0, Q“(s) > 0,0 < s < . 


lim P(- 0 n(r^ — r^) > x) = 1 — $ 
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And for any t < t, 


'^{'Tn <t) = ¥ (^Qn{s) < 0 or Qn{s) < 0, for some 0 < s < , 

= P (^Qnis) < 0 or Qn{s) < 0, for some 0 < s < . 

Now recall the large deviation principle for P(Q^(-) S £ •)> i-®-) Theorem fT^ and recall 

that f^{t) = + (j)^(t) — (j)‘^(t) — (j)^(t), Therefore, we have the 

following, 

Corollary 30. Given Theorem 1 1 91 for any t > t, 

lim — logP(r,i > t) = — inf I{4’)- 

^ g'^+0^(s)—0^(s)—0^(s)>O, 

for any 0<s<t 
0G^Co[O,oo) 

Similarly, for any t < t, 

lim — logP(r„ < t) = — inf ’^{4’)- 

n^oo n qb Jqj, gQjjig 0 < s < t 

or g°‘+4>'^{s)—<j>^(s)—(l>^{s)<0 for some 0 < s < t 
</ie.4Co[o,cxD) 


5 Extensions and discussions 


5.1 General assumptions for cancellation 

In the previous section, we have derived the fluid limit and fluctuations for the order positions 
under the simple assumption that cancellation is uniform on the queue. This assumption can be 
easily relaxed and the analysis can be modified fairly easily. 

For instance, one may assume (more realistically) that the closer the order to the queue head, 
the less likely it is cancelled. More generally, one may replace the term ( in Eqn. ( 12 . 21 ) with 

T where T is a Lipschitz-continuous increasing function from [0,1] to [0,1] with T(0) = 0 

and T(l) = 1. Now, the dynamics of the scaled processes are described as 


(Qi{t)\ /1 -1 

QUt) I = I 0 0 

\Zn{t) 


-1 0 0 0 
0 1 - 1-1 Iw 


0 -1 0 0 0 




■ddnit). 


(5.1) 


Then the limit processes would follow 

'Q\t) 

d I Q-{t) I = 


Z{t) 


/ 1 -1 -1 0 0 0 

0 0 0 1 -1 -1 

\o -1 -^( 1 !^) » 0 » 


^Q“(t-)>0,Q*(i-)>0,Z(t-)>0 ■ (V (b.2) 
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Theorem 31. Given Assumptions\^ [H and\^ and the scaled processes {QlQ-,Zr,) defined by 
Eqn. dSI]). If there exist constants q^, and z such that 

{Qi{0),Ql{0),Znm^{q\qfiz), 

then for any T > 0, Eqn. UM 

(Ql Q“, Z„) ^ (Q^ Q“, Z) in {D^[0, T], Ji), 
where (Q',Q“,z) is defined by Eqn. (1121) and 

(Q'>(0),Q“(0),Z(0)) = (g^(?^^). (5.3) 

Proof. First, let us extend the definition of T from [0,1] to M by 


T(x) = T(x)Io<a:<l + Il<x- 

Then T is (still) Lipschitz-continuous and increasing on M. That is, there exists K > 0, such 
that for any zi,Z 2 G K, |T(^;i) — T( 2 : 2 )| < K\zi — Z 2 \. Next, define r = min{r^,r“,r^} with 
= inf{t I Q^{t) < 0}, r“ = inf{t | Q°'{t) < 0}, and = inf{t | Z{t) < 0}. Similar to the 
argument for Lemma[ 6 l T e [0,1] and z,q’^ > 0 imply that Zn{t) < Q^(t) and Z{t) < Q^{t) for any 
time before hitting zero. Thus < r^. Now the remaining part of the proof is similar to that of 
Theorem [TT] except for the global existence and local uniqueness of the solution to Eqn. (15.2p . with 

^ = + (5.4) 

Denote the right hand side of Eqn. (15.4|) by '!?(Z, t), and define 'd{Z,q^/{Xv^)) = 1. Let {Ti}i>i be 
an increasing positive sequence with limj^oo= t- Then for any zi,Z 2 > 0 and 0 < t < Tj, 


|?9(zi,t) - '&{z2,t)\ = 


T 


< XV^K 


zi 

— Xv’^t^ 
zi 


-T 


^2 


^^2 


^q^ — Xv^t^ 


< 


XV^K 


q^ — Xv^t q^ — Xv^t 




Therefore '!?(Z, t) is Lipschitz-continuous in Z and continuous in t for any t <Ti and Z > 0. By the 
Picard’s existence theorem, there exists a unique solution to Eqn. (|5.4I) with the initial condition 
Z(0) = z on [0,Tj]. Now letting i —)• oo, the unique solution exists in [0,r). Moreover, by the 
boundedness of '&{Z,t) and the continuity of Z{t) at r, the unique solution also exists at t = r. 
For t > T, ^{Z, 0) = 0 and Z{t) = Z{t). Hence there exists a unique solution Z{t) for t > 0. Note 
that r“ = oo (resp. = oo) when < 0 (resp. < 0). However, since the right hand side of 
Eqn. ()5.4p is less than or equal to — it follows that Z{t) is decreasing in t and hits 0 in finite 
time. Therefore r is well defined. □ 
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5.2 Linear dependence between the order arrival and the trading volnme 

One may also replace Assumptions [T] and [2] by the assumption that order arrival rate is linearly 
correlated with trading volumes. The fluid limit can be analyzed in a similar way with few modi¬ 
fications. 

Assumption 32. N{nt) is a simple point process with an intensity nX + anQ'^ (t— ) -)- (3nQ^ (t— ) 
at time t, where a, /3 are positive constants. 

Assumption 33. For any 1 < j <Q, is a sequence of stationary, ergodic, and uniformly 

bounded sequence. Moreover, for any i>2 and Qi = <7{{^ k}i<k<i), 

I = y. 


Theorem 34. Given Assumptions El EH and EH then Theorem [17] holds except that the limit 
processes will be replaced by 


Q\t) = 


aq^v^ — aq^v°' + \v^ 
v°-a -\- 


^ v'^{l3q^ + aq°- + A) ^-(vb 0 j^yo.a)tAT 
jdv^ + av°- 


(5.5) 


^ - I3q''v + v'^{/3q^ + aq^ + A) (^bg+^aa)tAT 

v°-a -|- v^j5 (3v^ + av°- ’ 


and 


Z{t) = ze 




Q»{s) ' ^ • Q>>(s) 


ds 


ptAr _ _ t 7 r A I 4 ? I c>:Q^(u) ~| 

- / V 2 [X + l3Q\s) + aQ'^{s)]e ds. 

Jo 

Proof. Recall that before t < t, with Assumption 1321 


(Qiit)' 

crjt) 1 = 10 0 

\Zn(*) 


1-1 -1 000 
0 1 -1 -1 11 , 

“ -llwj » » “ 


Q“(t-)>0,Z„(t-)>0 


• d^nit), 


where 


Here 


^ Mint) .t 

'A n{t) = — 'y ^ V i = Mn{t) + / {X + PQnis — ) + 0'Q^{s — ))dsV. 
” i=l -^0 


'^ln(t) = 


N{nt) 


(5.6) 


(5.7) 


W = - + N{nt)-n[ (A + /3Q^(s-) + aQ“(s-))d 

n ^ n [ Jo 

is a martingale. Similar to the arguments before, we can show that (QtQ“,Z„) ^ (Q',Q“,Z) 
where (Q^,Q“,Z) satisfies the ODE: 


(Q\ty 

QAt) I = 
\z{t) 


/ 1 -1 -1 0 0 0 




0 0 0 1 -1 -1 I + + 

0 -1 0 0 0 
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with the initial condition ((5^(0), (5“(0), Z(0)) = {q’^,q°',z). The equations for Q’^{t) and Q°'{t) can 
be written down more explicitly as 


dQ\t) = (A + /3Q\t-) + aQ<^{t-)){Vi -V 2 - Vz)dt, 
dQ^{t) = (A + PQ\t-) + aQ^{t-)){Vi -V^- VG)dt, 


which can be further simplified as 


d 


Q\t) 

Q^{t) 


Hence, for t < r, we get 


-u^/3 -v^a \ ( Q\t) \ 
—u“/3 —v^a J y Q“(t) / 


( Q\t) 

V 


Cl 


a 

-/3 


+ C2e 


— {v^P+v°'a)t 


Xv^ \ 
Au“ J ’ 



where ci, C2 are constants that can be determined from the initial condition, 

^ v^a + u^/3 ’ ^ [5v^ + av°- 

Hence Eqns ()5.5I) and (15.61) follow. 

Finally, Z{t) satisfies the first-order ODE 

dZ{t) + Z{t)Vs + /3 + ^11^) dt = -V2{X + fiQ\t) + aQ-{t))dt, 

whose solution is given by Eqn. dSZD- □ 

Corollary 35. Given Assumptions\^ EH and EH assume further that v^/3 + u“a > 0 and < 
qayb _ gbya ^ Then Q^{t) and (5“(t) will hit zero at some finite times and r“ respectively. 
Moreover, 


J 1 f v^X + q°‘v^a — q^v^-aX 

+ v^-a ^ \ v^fiq^ + v^aq^ + Xv^ J ’ 

^ 1 /—q°‘v'^/3 + + Xv°‘\ 

+ v°‘a ^ y -|- aq^-v^- + Xv°- J ’ 

and is determined via the equation 


z = 


f 


H2(A + l5Q\s) + aQ“(s))e^° J ds 


foVs 


T+/3+ 




du 


5.3 Various forms of diffusion limits 

There is more than one possible alternative set of assumptions under which appropriate forms of 
diffusion limits may be derived. For instance, one may impose a weaker condition than Assumption 
[12] for {A}i>i- 
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Assumption 36. For any time t, 

Nint) 

lim -= At, a.s. 

n—^oo Tl 

Moreover, there exists K > 0, such that E[A(t)] < Kt, for any t. 

This assumption holds, for example, if the point process N{t) is stationary and ergodic with 
finite mean. To compensate for the weakened Assumption 1361 one may need a stronger condition 
on for instance, Assumption 1331 

Note that under this alternative set of assumptions, the resulting limit process will in fact be 
simpler than Theorem 1 141 This is because Assumption 1331 implies that V/ is actually uncorrelated 
to V^, for any i ^ i' and 1 < j < 6. Hence the covariance of and V/, i > 2 m the limit process 
may vanish. We illustrate this in some detail below. 

Making Assumptions [33] and (361 define a modified version of the scaled net order flow process 

by 

N{nt) 

¥;(t) = — ^ (V^-V), (5.8) 

while the scaled processes Rn{t), Rfit) still follows Eqn. (13.3p . the hrst hitting time the same as 
in Eqn. (13.5p . and the corresponding limit processes in Eqn. ()3.13p and Eqn. (I3.12p . Then we have 
the following. 

Theorem 37. Given Assumptions\^ \3^ and\3^ for any T > 0, 

(i) ^ where = {ajWj,l < j < 6), where {Wj,l < j < 6) is a standard six¬ 

dimensional Brownian motion and Uj = XVariyl). 

(R^,R“)^(R',R“) in 

Proof. Under Assumption 1331 it is clear that 

N{nt) 

V" i=i 

is a martingale. Now define for j = 1, 2,..., 6, 

N{nt) N{nt) 

ML ■■= E iy' - I = E (y - 

Eirst, the jump size of is uniformly bounded since N{nt) is a simple point process and by 
Assumption [33l VFs are uniformly bounded. Next, the quadratic variation of is given by 

W (nt) 

2=1 

By Assumptions [331 and 1361 and the Ergodic theorem, as t —)• oo, 

—)• AVar]!/-^], a.s. 
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Moreover, since JVP and have no common jumps for j ^ k, 


= 0 . 

Therefore, applying the FCLT for martingales m Theorem VIII-3.11], for any T > 0, we have 

m{D%T],J,), 

To see the second part of the claim, first note that by Assumption 1361 


Af(n-) ^ ^ 

-Vy^AFe, m(D[0,T],Ji) a.s. 
n 

1 = 1 


as n —oo. The remaining of the proof is to check the conditions for Theorem 1101 as in the proof of 
Theorem m The quadratic variance of Mnt '■= is given by 


E 


\-^M 



nt- 


n 


E[iV(nt)]E 


i<i<6 


Vj - E 


vj I j; 


i I T/- 


<Kt Y ^ 

i<i<6 




which is uniformly bounded in n. The total variation of An := ^ ^ satisfies 




n[TiAn)]t]< Y^ 

i<i<6 


N(nt) 


-n E 


2=1 


< KtE[\V^\], 

i<i<6 


which is uniformly bounded in n. The proof is complete. 


□ 


References 

[1] F. Abergel and A. Jedidi. A mathematical approach to order book modeling. International 
Journal of Theoretical and Applied Finance, 16(05), 2013. 

[2] F. Abergel and A. Jedidi. Long time behaviour of a Hawkes process-based limit order book. 
Preprint, 2015. 

[3] A. Alfonsi, A. Fruth, and A. Schied. Optimal execution strategies in limit order books with 
general shape functions. Quantitative Finance, 10(2): 143-157, 2010. 

[4] A. Alfonsi, A. Schied, and A. Slynko. Order book resilience, price manipulation, and the 
positive portfolio problem. SIAM Journal on Financial Mathematics, 3(l):511-533, 2012. 

[5] S. Asmussen and H. Albrecher. Ruin Probabilities. World Scientific, Singapore, 2010. 

[6] M. Avellaneda and S. Stoikov. High-frequency trading in a limit order book. Quantitative 
Finance, 8(3):217-224, 2008. 

[7] E. Bayraktar and M. Ludkovski. Liquidation in limit order books with controlled intensity. 
Mathematical Finance, 24(4):627-650, 2014. 


35 














[8] P. Billingsley. Convergence of Probability Measures. John Wiley, New York, 1968. 

[9] J. Blanche! and X. Chen. Continuous-time modeling of bid-ask spread and price dynamics in 
limit order books. arXiv preprint arXiv:1310.1103, 2013. 

[10] R. C. Bradley. Basic properties of strong mixing conditions, a survey and some open questions. 
Probability Surveys, 2:107-144, 2005. 

[11] L. Breiman. Probability. Classics in Applied Mathematics. Society for Industrial and Applied 
Mathematics, 1968. 

[12] P. Bremaud and L. Massoulie. Stability of nonlinear Hawkes processes. The Annals of Proba¬ 
bility, 24(3): 1563-1588, 1996. 

[13] W. Bryc. On large deviations for uniformly strong mixing sequences. Stochastic Processes and 
their Applications, 41 (2): 191-202, 1992. 

[14] W. Bryc and A. Dembo. Large deviations and strong mixing. Annales de I’lHP Probabilites 
et statistiques, 32(4):549-569, 1996. 

[15] A. Bulinski and A. Shashkin. Limit theorems for associated random fields and related systems. 
World Scientific, Singapore, 2007. 

[16] R. M. Burton, A. Dabrowski, and H. Dehling. An invariance principle for weakly associated 
random vectors. Stochastic Processes and their Applications, 23(2):301-306, 1986. 

[17] A. Cartea and S. Jaimungal. Modelling asset prices for algorithmic and high-frequency trading. 
Applied Mathematical Finance, 20(6):512-547, 2013. 

[18] A. Cartea, S. Jaimungal, and J. Ricci. Buy low, sell high: A high frequency trading perspective. 
SIAM Journal on Financial Mathematics, 5(l):415-444, 2014. 

[19] R. Cont and A. De Larrard. Order book dynamics in liquid markets: limit theorems and 
diffusion approximations. Available at SSRN 1757861, 2012. 

[20] R. Cont and A. De Larrard. Price dynamics in a Markovian limit order market. SIAM Journal 
on Financial Mathematics, 4(1): 1-25, 2013. 

[21] R. Cont and A. Kukanov. Optimal order placement in limit order markets. Available at SSRN 
2155218, 2013. 

[22] R. Cont, S. Stoikov, and R. Talreja. A stochastic model for order book dynamics. Operations 
research, 58(3):549-563, 2010. 

[23] A. Dembo and T. Zajic. Large deviations: from empirical mean and measure to partial sums 
process. Stochastic Processes and their Applications, 57(2): 191-224, 1995. 

[24] A. Dembo and O. Zeitouni. Large Deviations Techniques and Applications. Springer, New 
York, 1998. 

[25] P. W. Glynn and W. Whitt. Ordinary CLT and WLLN versions of L=AW. Mathematics of 
operations research, 13(4):674-692, 1988. 


36 


[26] O. Gueant, C.-A. Lehalle, and J. Fernandez-Tapia. Optimal portfolio liquidation with limit 
orders. SIAM Journal on Financial Mathematics, 3(l):740-764, 2012. 

[27] F. Guilbaud and H. Pham. Optimal high-frequency trading with limit and market orders. 
Quantitative Finance, 13(l):79-94, 2013. 

[28] X. Guo. Optimal placement in a limit order book. TUTORIALS in Operations Research, 
INFORMS, 2013. 

[29] X. Guo, A. De Larrard, and Z. Ruan. Optimal placement in a limit order book: an analytical 
approach. Preprint, 2013. 

[30] A. G. Hawkes. Spectra of some self-exciting and mutually exciting point processes. Biometrika, 
58(l):83-90, 1971. 

[31] U. Horst and D. Kreher. A weak law of large numbers for a limit order book model with fully 
state dependent order dynamics. arXiv preprint arXiv: 1502.04359, 2015. 

[32] U. Horst and M. Paulsen. A law of large numbers for limit order books. arXiv preprint 
arXiv:1501.00843, 2015. 

[33] W. Huang, G.-A. L. Lehalle, and M. Rosenbaum. Simulating and analyzing order book data: 
The queue-reactive model. Journal of the American Statistical Association, 110 (509):107~122, 
2015. 

[34] H. Hull and J. Kiessling. Algorithmic trading with Markov chains. PhD thesis. Doctoral 
Thesis, Stockholm University, Sweden, 2010. 

[35] I. A. Ibragimov. A note on the central limit theorem for dependent random variables. Theory 
Prohah. Appl, 20:135-140, 1975. 

[36] S. Iyengar. Hitting lines with two-dimensional Brownian motion. SIAM Journal on Applied 
Mathematics, 45(6):983-989, 1985. 

[37] J. Jacod and A. N. Shiryaev. Limit Theorems for Stochastic Processes. Springer Berlin, 1987. 

[38] A. Kirilenko, R. B. Sowers, and X. Meng. A multiscale model of high-frequency trading. 
Algorithmic Finance, 2(l):59-98, 2013. 

[39] L. Kruk. Functional limit theorems for a simple auction. Mathematics of Operations Research, 
28(4):716-751, 2003. 

[40] T. Kurtz and P. Protter. Weak limit theorems for stochastic integrals and stochastic differential 
equations. The Annals of Probability, 19(3):1035-1070, 1991. 

[41] S. Laruelle, G.-A. Lehalle, and G. Pages. Optimal split of orders across liquidity pools: a 
stochastic algorithm approach. SIAM Journal on Financial Mathematics, 2(1): 1042-1076, 
2011. 

[42] G. Maglaras, G. G. Moallemi, and H. Zheng. Optimal order routing in a fragmented market. 
Preprint, 2012. 


37 


[43] A. Metzler. On the first passage problem for correlated Brownian motion. Statistics & proba¬ 
bility letters, 80(5):277-284, 2010. 

[44] C. C. Moallemi and K. Yuan. The value of queue poistion in a limit order book. Working 
paper, 2015. 

[45] S. Predoiu, G. Shaikhet, and S. Shreve. Optimal execution in a general one-sided limit-order 
book. SIAM Journal on Financial Mathematics, 2(1):183-212, 2011. 

[46] M. Rosenblatt. Markov Processes, Structure and Asymptotic Behavior. Springer-Verlag, New 
York, 1971. 

[47] S. E. Shreve, C. Almost, and J. Lehoczky. Diffusion scaling of a limit-order book model. 
Working paper, 2014. 

[48] C. Tone. A central limit theorem for multivariate strongly mixing random fields. Probab. 
Math. Statist, 30(2):215-222, 2010. 

[49] S. R. S. Varadhan. Large Deviations and Applications. SIAM, Philadelphia, 1984. 

[50] L. A. Veraart. Optimal market making in the foreign exchange market. Applied Mathematical 
Finance, 17(4):359-372, 2010. 

[51] A. R. Ward and P. W. Glynn. A diffusion approximation for a Markovian queue with reneging. 
Queueing Systems, 43(1-2):103-128, 2003. 

[52] A. R. Ward and P. W. Glynn. A diffusion approximation for a GI/GI/1 queue with balking 
or reneging. Queueing Systems, 50(4):371-400, 2005. 

[53] W. Whitt. Stochastic-Process Limits: An Introduction to Stochastic-Process Limits and their 
Application to Queues. Springer, New York, 2002. 

[54] G. Zhou. An analysis of default correlations and multiple defaults. Review of Financial Studies, 
14(2):555-576, 2001. 

[55] L. Zhu. Gentral limit theorem for nonlinear Hawkes processes. Journal of Applied Probability, 
50(3):760-771, 2013. 

A Some large deviations results 

According to |23l Theorem 2], we have 

Theorem 38. Let be a sequence of stationary -valued random vectors satisfying As- 

sumption{T^ and AssumptionlT7[ Then, the empirical mean process Sn{t) := ^ Yli=i 0 < t < T, 

satisfies a large deviations principle on D[0,T] equipped with the topology of uniform convergence 

with the convex good rate function 



(A.l) 
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for any cp £ -4.Co[0, oo), the space of absolutely continuous functions starting at 0 and T{4>) = +oo 
otherwise, where 

K{x) ■.= sup {0 • X — r(0)}, (A.2) 

with T{9) := lim^^oo ^ logE[e^^=i®'^*]. 

Remark. Note that the original statement in [231 Theorem 2] applies to Banach space-valued 
For the purpose in our paper, we only need to consider R^-valued 


B Y(t) process 

Proposition 39. Y{t) defined in Eqn I13.16\} is a Gaussian process for t < , with mean 0 and 

variance (7y(t). In particular, when c < 0 and c ^ —1, 
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with 


a := -('012 - 022 - 032) + (013 - 023 “ 033 ) 
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Remark. Proposition [39] only gives the formula for the variance of Y{t) for the case c 7 ^ — 1, 
c < 0. The variance cTy(t) for the case c = —1 can be taken as a continuum limit as c ^ —1. 

Proof of Proposition (23 By multiplying Eqn. (I3.16p by the integrating factor e'^° and inte¬ 

grating from 0 to t, and finally dividing the integrating factor, we get 


rt Av^ j r 

Y{t) = - e ^ - 

Jo Jo 

I 


0 Q\s) 


(B.3) 


* -S: ^SuZ{s)i^\s) - ^Hs) - T3(s)) ,-3 


+ / e *3 


XV-^ds, 


which implies that Y(t) is a Gaussian process since ^ is a Gaussian process. Since ^ is centered, 
i.e., with mean zero, it is easy to see that Y{t) is also centered. Next, let us determine the variance 
of Y{t). By Ito’s formula, we have 


d{Yitf) = 2Y{t)dY{t) + d{Y)t 

= d{Y)t - 2Y{t)^^XV^dt - 2Y{t)d^\t) - 2y(t)^dT3(t) 


(B.4) 


From Eqn. (|3.16l) . we get 


d{Y)t = d{^^)t + 


2\ I j/,u3\ , 2Z(t) 


QHt? 




(B.5) 


Plugging Eqn. (IB.Sp into Eqn. (IB.4p . and taking expectations on the both hand sides of the equation, 
we get 

dE[Y{tf] = d{^% + -^d{^^)t + 
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rt 2XV^ 

By using the integrating factor e'^° , we conclude that 
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Let us recall that 
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We also recall that (V'ij)i<i,j<6 is a symmetric matrix defined as 


V’ij := ^ + V-V^vjX^, l<i,j <6. 


k=l 


Therefore, we have 


For any i,j and t > s, 


= '4^22t, (^^)t = V'sai, (^^ = ^>23^ 


E[T'(t)T^(s)] = ^ S.fcSjfcAs + V^V^vlXh = 
k=l 

For any i = 1,2,3, from Eqn. (IB.3p . we can compute E[y(t)'I'*(t)] as 
E[y(t)T*(t)] 
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Next, combining Eqns. (IB.7p . (IB.Sp . (IB.9h . (I2.17p . and (I2.1ip . and after some calculations, we 
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where a,/3,j,S are defined in Eqn. (IB.2I) and are defined in Eqn. (IB.II) . Therefore, 
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(B.12) 


Hence, we get the desired result by substituting Eqn. (IB.lOp and Eqn. (IB.1211 into Eqn. (IB.6E □ 
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